Electronic structure of epitaxial graphene layers on SiC: effect of the substrate. 
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Recent transport measurements on thin graphite films grown on SiC show large coherence lengths 
and anomalous integer quantum Hall effects expected for isolated graphene sheets. This is the 
case eventhough the layer-substrate epitaxy of these films implies a strong interface bond that 
should induce perturbations in the graphene electronic structure. Our DFT calculations confirm 
this strong substrate-graphite bond in the first adsorbed carbon layer that prevents any graphitic 
electronic properties for this layer. However, the graphitic nature of the film is recovered by the 
second and third absorbed layers. This effect is seen in both the (OOOl)and (0001) 4H SiC surfaces. 
We also present evidence of a charge transfer that depends on the interface geometry. It causes the 
graphene to be doped and gives rise to a gap opening at the Dirac point after 3 carbon layers are 
deposited in agreement with recent ARPES experiments (T.Ohta et al, Science 313 (2006) 951). 

PACS numbers: 73.20.At, 71.15.Mb 
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The possibility of carbon nanotubes (CNT) switching 
devices has been pursued in the last decade because of 
their attractive electronic properties. Nevertheless, prob- 
lems with large intrinsic resistance in contacts and the 
inability to control tube helicity, and thus whether or 
not they are metallic or semiconducting, have made large 
scale circuit designs problematic. The proposed solu- 
tion to these problems is an all carbon nanoelectron- 
ics paradigm based on the planar 2D form of carbon, 
graphene. [l| 

Graphene consists of a single carbon plane arranged on 
a honeycomb lattice. From a fundamental point of view, 
graphene ribbons can be seen as an unrolled CNT but 
with different boundary conditions (finite versus cyclic). 
Therefore, their electronic properties should be similar. 
In fact this has been demonstrated in recent experiments 
on single and multi-graphene sheets that show the exis- 
tence of Dirac Fermions, large electron coherence lengths 
and anomalous integer quantum Hall effect 0, 0, 0). 
The advantage of graphene over CNTs for electronics re- 
sides in its planar 2D structure that enables circuit de- 
sign with standard lithography techniques. This enables 
the graphene to be cut with different shapes and selected 
edge direction. By simply selecting the ribbon edge di- 
rection it is possible to design metallic or semiconductor 
graphene ribbons (analogous to helicity in CNTs). 

Since single or multiple sheets must be supported on 
a surface for fabrication, the pressing question becomes: 
how does the interface between a graphene sheet and its 
support affect its electronic properties? In other words 
can the symmetry of an isolated graphene sheet be main- 
tained in the presence of an interface? It is this question 
that is the focus of this paper. Specifically we have stud- 
ied the system of graphite grown on both polar faces of 
hexagonal SiC. 

The graphene layers are produced by sublimating Si 



from either the 4H- or 6H-SiC (0001) (Si terminated) or 
(0001) (C terminated) surfaces at sufficiently high tem- 
peratures to graphitize the excess carbon [jj, [7| . 2D trans- 
port measurements on these graphitized surfaces show 
the presence of Dirac electrons similar to those found on 
exfoliated graphene [2, H, H, @| . Besides being a more 
practical and scalable approach to 2D graphene electron- 
ics, this system has the experimental advantage of having 
a well defined interface that can be characterized in con- 
trast to mechanically exfoliated graphene flakes [H, 0| 
that must in any case still be supported on a surface 
(usually SiC^)- 

In this letter, we conclusively show that the first car- 
bon layer grown above the SiC substrate has no graphitic 
electronic properties and acts as a buffer layer between 
the substrate and subsequent graphene layers. Atoms in 
this plane form strong covalent bonds with the SiC sub- 
strate. The next graphene layer above the buffer layer 
shows a graphene-like Dirac band structure expected for 
an isolated graphene sheet. The calculated results are 
consistent with a short C first plane-4HSiC(0001) sub- 
strate bond as determined by Surface X-ray Reflectivity. 
Under some conditions, charge transfer from the sub- 
strate results in a n type doping of the graphene layers 
(Fermi level above the Dirac point). This opens a gap 
in the graphene bilaycr Dirac bands, in agreement with 
recent ARPES results [l0( . Dangling bond related states 
are found in all tested geometries. These states can in- 
teract with the graphene derived state depending on the 
geometry. Their effects on the electronic and transport 
properties have to be considered as well as those of the 
intrinsic defects of the isolated graphene layer [llj . 

The systems theoretically studied here are made of one, 
two or three carbon layers (on a honeycomb lattice with 
Bernal stacking) on top of either a SiC 4H (0001) or 
(0001) (Si and C terminated) substrates. Graphene is 
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nearly commensurate with these SiC surfaces with a com- 
mon cell corresponding to a 6-\/3x6v / 3~R30 (with respect 
to SiC lxl surface cell) @- This cell is too large to 
make realistic calculations. Even the next smallest near- 
est commensurate structure, 4x4 cell (not experimentally 
observed), is too large for reasonable calculations. There- 
fore, as a first approximation to the actual structure, we 
used the \/3x\/3-R30 reconstruction shown in Fig.Q] The 
y3x V3R30 cell corresponds to a 2x2 graphene cell. The 
graphene and SiC lattice parameter mismatch requires an 
8% stretch of the graphene to make the two cells com- 
mensurate. We have checked that this expansion has no 
qualitative effect on a free standing graphene electronic 
structure (it will however change the Dirac electron ve- 
locity). The cell contains 3 atoms/layer in SiC. At the 
interface, two of these atoms are immediately below a C 
atom in the first C layer. The third atom (subsequently 
referred to as the "lonely atom") has no C atom above 
it. A bulk truncated 4H-SiC geometry was used on both 
faces. We have also checked another possible interface 
geometry based on surface X-ray scattering data: a C- 
terminated surface with one C atom out of three missing 
(referred to as " C-deficient" ) . In the bulk truncated ge- 
ometry, the lonely atom exhibits a dangling bond (DB) 
that points towards the graphene layers. The lonely atom 
is suppressed in the C-deficient geometry, thus creating 
3 dangling bonds at the interface. 

The electronic structure was investigated using the 
VASP code fl2|| . It is based on Density functional the- 
ory (DFT) within the generalized gradient approximation 
[13J . The 4H SiC substrate is described with a slab that 
contains 8 SiC bi-layers with H saturated dangling bonds 
on the second surface. The empty s pac e ranges from 15 
to 25 A. Ultra soft pseudopotentials 14 1 are used with a 
plane wave basis cutoff equal to 211 eV. The experimen- 
tal graphene layer spacing was first chosen as the starting 
value and then all the atoms were allowed to relax. Since 
DFT is known to poorly describe Van der Waals forces, 
the final graphene layer spacing are significantly larger 
than in the bulk. However, we point out that the C-short 
ultrasoft pseudopotential used here has been extensively 
tested [15j and was shown to correctly describe the band 
structure of graphite in spite of the larger layer spacings 
lil T3] • Integration over the Brillouin zone is performed 
on a 9x9x1 grid in the Monckhorst-Pack scheme to en- 
sure convergence of the Kohn Sham eigenvalues. The K 
point is included since it is crucial to a good description 
of the Fermi level for a single graphene layer. 

The X-ray experiments were performed at the Ad- 
vanced Photon Source, Argonne National Laboratory, on 
the 6IDC-/iCAT Ultra High Vacuum scattering chamber. 
The C-Face samples were graphitized in a vacuum RF- 
induction furnace (P — 3 x 10~ 5 Torr) [l8[ and transported 
to the scattering chamber for analysis. 

The calculations show that the relaxed geometry of the 
bulk layers is influenced only by the first carbon layer. 



Neither the bulk nor the first C layer are altered when 
subsequent C layers are added. In the last bilayer of 
the Si (C) terminated face the lonely atom is displaced 
toward the bulk by 0.3 A (0.45), while the bilayer width 
remains globally unchanged (0.65 (0.7) compared to 0.62- 
0.63 in the bulk). The first graphene layer lies 2.0 (1.66) 
A above the two outermost atoms. The second graphitic 
C plane is 3.8 A (3.9) above the first one (subsequent 
planes are spaced by 3.9 A). As mentionned above, the 
large value of the graphene interlayer spacings is due to 
the known difficulty of the DFT to describe the van der 
Waals force [lj| 17 1. From these results, we deduce that 
the interface carbon layer strongly interacts with the SiC 
substrate for both Si and C terminated surfaces. Sub- 
sequent C-planes on the other hand are weakly bound 
by van der Waals forces as expected for graphite. This 
conclusion also holds for the C-deficient geometry. The 
bulk relaxation in this latter case is very similar to those 
of the C terminated bulk truncated geometry. 

X-ray reflectivity data confirms this result. Figure. [2] 
shows an experimental reflectivity from ~ 9 graphene 
layer film grown on the 4H-SiC(0001) surface. Data is 
on the (00Z) rod in units of 27r/a, where a = 10.08lA. 
Peaks at I = 4 and 8 are SiC Bragg reflections, while 
peaks at I ~3, 6 and 9 are graphite Bragg points. A 
full fit to the data including substrate relaxation and a 
multilayered graphite film is shown. Details of the fit 
are given in a separate article The model consists 
of a single reconstructed SiC bi-layer interface between 
the graphite and the bulk. The fit reveals that the first 
graphene layer is 1.65 ± .05A above the last bulk C-layer 
consistent with the calculated value. The next graphene 
layer is separated from the first by 3.51 ± O.lA (slightly 
larger than the bulk value 3.354A [2(|). Subsequent lay- 
ers have a mean spacing of 3.370 ± 0.005A. This slightly 
larger layer spacing is consistent with stacking faults in 
the layers [2lJ. The x-ray results confirm the calculated 
structure of a strongly bonded first graphitic layer with 
a well isolated graphene layer above it. We note that the 
extended diamond interface phase conjectured by others 
0, IS] does not fit the x-ray data for the C-face. This 
is demonstrated in Fig. [2] were we force a second SiC bi- 
layer to be Si depleted by 25%. This fit is obviously worse 
than a single bi-layer and proves that the SiC interface 
is narrow and not extended. 

The band structures with one (Fig.[3ja) and (b)), two 
(Fig. EJc) and (d)) and three (Fig. [3](e) and (f)) carbon 
layers on bulk terminated SiC are shown in Fig. [3] For 
both polarities, the electronic structure with a single C 
layer significantly differs from graphite [23[. It exhibits 
a large gap and a Fermi level pinned by a state with a 
small dispersion (close to the conduction band for the Si 
terminated surface or in the gap for the C terminated 
surface). These states are related to the Dangling Bond, 
DB, of the lonely atom in the SiC interface layer (a Si (C) 
dangling bond state for the Si (C) terminated surface). 
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They remain unchanged when further C layers are added 
on top of the first one. Figures. [3jc-f) show that graphene 
related dispersions are recovered when more than one 
C layer is present. In fact, the first C plane acts as a 
buffer layer that allows growth of subsequent graphene 
like layers. Indeed one can clearly see for two C planes 
(buffcr+1) the linear dispersion and Dirac point that are 
characteristic of an isolated graphene layer (Fig. [3]^c-d)). 
When three C layers are present (buffcr+2) the disper- 
sion is similar to the dispersion of a graphene bilayer 

On Si-terminated surface, the Fermi level falls 0.4 eV 
above the Dirac point. The graphene like planes are n 
doped. This is confirmed by the opening of a gap in the 
case of 3 C layers (buffer +2). Tight binding calculations 
involving p z orbitals show that this is characteristic of a 
graphene bilayer where the two planes are not symmetric. 
In our ab initio calculation, the Fermi level falls above the 
highest unoccupied 7r* band minimum at K point. The 
comparison to tight binding calculations, shows that the 
two graphene layers are doped and that one plane is less 
doped than the other one. This is in agreement with 
recent ARPES and XPS measurements jl, 03 • For C 
terminated surfaces in the bulk truncated geometry, the 
Fermi level falls on the Dirac point and the graphene 
layers are neutral. On the other hand, the Fermi level 
of the C-deficient geometry (Fig. is 0.4 eV above the 
Dirac point. It is fixed by states related to the 3 DB 
present in this structure. This stresses the role played by 
interface defects. 

For the Si terminated surface the DB related state and 
the graphene derived states anticross indicating some in- 
teractions between them (inset of Fig. 3c). This is not 
the case for the bulk C terminated surface. This effect 
may have a crucial impact on transport properties of the 
film and explain the low electron mobilities of Si-face 
films compared to C-face films. f5~s(] 

In Fig. [5j charge density contours show clear evidence 
of the existence of a covalent bond between the first 
graphitic C layer (buffer layer) and SiC. Charge den- 
sity appears to be more delocalizcd in the subsequent 
graphitic planes as one can expect for graphene like lay- 
ers. 

In conclusion, we have shown that the first C layer 
on top of a SiC surface acts as a buffer layer and allows 
the next graphene layer to behave electronically like an 
isolated graphene sheet. The existence of strong cova- 
lent bonds to the first layer is in agreement with X-ray 
reflectivity data. The electronic structure of subsequent 
graphitic layers depends on the geometry of the interface 
and on the number of layers. Our calculations show clear 
evidence of a charge transfer from SiC to the graphene 
layers that depends on the interface geometry and results 
in a doping of these layers. We even show the possible 
opening of a gap at the Dirac point in agreement with 
ARPES results. Interface intrinsic defects induce states 



in the vicinity of the Fermi level. The interaction of these 
states with the graphene derived states depends on in- 
terface geometry and may explain the lower electronic 
mobility observed on Si-terminated surface. Because the 
defect density (i.e. DB states) is even larger for the actual 
6-^/3 x6-\/3-R30 cell, further experiments and calculations 
are needed to clarify the role of these states and their 
dependence on interface geometry and stacking order in 
these systems. 
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FIG. 2: (color on line)X-ray (00Z) reflectivity data from 9 
graphene layers grown on the 4H-SiC(0001) surface. Bulk 
and graphite Bragg peaks are labelled. Blue line is the best 
fit structure with one reconstructed SiC bi-layer as described 
in the text. Red line is a fit with an extended reconstruction 
of two SiC bi-layers. 



FIG. 1: (color on line) interface geometry, a- side view, b- top 
view of the \/3 x V3~R30 cell in the case of a Si-terminated SiC 
face. The lonely atom is missing in the C-deficient geometry 
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FIG. 3: (color on line) dispersion curves for one (a,b), two 
(c,d) and three (e,f) C layers on bulk truncated SiC. a,c,e 
correspond to Si terminated face; b,d,f to C terminated face. 
Inset in c shows a zoom of the anticrossing in the vicinity of 
E F . 




FIG. 4: (color on line) dispersion curves for 2 C layers on top 
of the C-deficient surface 




FIG. 5: (color on line) isocharge density contours along z axis 
for three C layers on C- (a) and Si-(b) terminated surface 



